#ifndef _EVAL_UTIL_H_
#define _EVAL_UTIL_H_

#include <surftrack.h>

using namespace eltopo2d;

void assemblyBoundaryVector(SurfTrack* surf, double** b);
void myMatrixVectorProduct(double* x, double** y, int n);
void evaluateDisplacementField(double* theta, double** u);
void evaluateDeformationGradientField(double* theta, double** F);


#endif